Transient Shear Banding in Time-dependent Fluids 
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We study the dynamics of shear-band formation and evolution using a simple rheological model. 
The description couples the local structure and viscosity to the applied shear stress. We consider in 
detail the Couette geometry, where the model is solved iteratively with the Navier-Stokes equation 
to obtain the time-evolution of the local velocity and viscosity fields. It is found that the underlying 
reason for dynamic effects is the non-homogeneous shear distribution, which is amplified due to a 
positive feedback between the flow field and the viscosity response of the shear thinning fluid. This 
offers a simple explanation for the recent observations of Transient Shear Banding in time-dependent 
fluids. Extensions to more complicated rheological systems are considered. 



The property that characterizes complex fluids is their 
non-trivial rheology, shear rate - stress relation. They 
are generally further categorized into shear thinning or 
shear thickening fluids. Both cases are additionally com- 
plicated by time dependence. Due to the stress-shear in- 
teraction, already small perturbations in the local stress 
can result in a positive feedback with the flow promot- 
ing shear instabilities in each case [H Q]- The under- 
standing of complex fluids is of enormous importance for 
many practical applications [3i] and the theory touches 
on many branches of physics. Recent advances allow to 
follow the suspension local velocity during a standard 
rheological experiment 0, @| ■ Quantifying the local flow 
field simultaneously with rheological measurements gives 
the possibility to measure both the intrinsic and apparent 
rheology. This has lead to the discovery that a hetero- 
geneous shear distribution in samples during such tests 
is ubiquitous. Shear banding [6] has been observed in 
many systems composed of substantially different build- 
ing blocks, such as colloidal glasses, wormlike micelles, 
foams and granular matter [7]. The current viewpoint, 
both phenomenologically and theoretically, is that a non- 
monotonic intrinsic flow curve, is what is common to 
most of these materials Q , but also other mechanisms 
have been suggested 

A branch of complex fluids are the simple yield stress 
fluids These materials do not show aging phenom- 
ena (thixotropy) . Therefore they are expected to have a 
monotonic intrinsic flow curve, and a steady state with- 
out shear bands [lH. However, recent experiments (l2| 
display shear banding during startup flows in a rotational 
rheonieter indicating time-dependent behavior. These, 
so called transient shear bands, can be very long last- 
ing, but eventually vanish with a homogeneous steady 
state. The transient shear banding phenomenon tests 
our fundamental understanding of non- Newtonian fluids, 
and is also important for industrial processes and sim- 
ply for understanding usual rheological measurements. A 
particular feature of the transient shear banding is that 
it appears to exhibit scaling familiar from critical phe- 



nomena: the time it takes for the transient to disappear 
(fluidization time r/) is a power-law function of the shear 
rate or applied stress [l^ . 

For the reasons leading to transient shear banding, see 
the summary by Adams et al. (l3| . three main candi- 
dates are offered: i) with certain parameters, the dy- 
namical equations are unstable amplifying small pertur- 
bations, which slowly quench towards the homogeneous 
steady state ii) The fluid flow curve is time depen- 
dent, and can, at different times, have non-monotonic 
shape [31 . iii) Elastic stress overshoots cause instabil- 
ity in the flow [l^. Some theoretical models appear to 
produce transient shear banding including shear trans- 
formation zone theories a modifled soft glassy rhe- 
ology model [l^, a simplifled fluidity model and a 
mesoscopic model of plasticity |19,]. Such models of tran- 
sient shear banding share the property of time-dependent 



17 



reduction of the local stress under shear 
plained in Ref. [2^. Literature reports experimental ev 
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idence of time-dependent rheology in carbopol gels [12 1 
especially at small shear rates, as well as in other simple 
yield stress fluids [2l|, appearing as slight hysteresis in 
the flow curves. Further details of the time-dependent 
flow curve hysteresis related to carbopol gels is reported 
in Ref. |22l| . Stronger hysteresis is inherent to thixotropic 
fluids [21|. Indeed, transient shear banding has been re- 
cently found also there well above the critical shear rates 
(2^ indicating that transient shear banding could be a 
general property of soft glassy materials. 

Motivated by these findings, we consider here a struc- 
tural model to find the main ingredients of transient shear 
banding. It is for a simple time-dependent Newtonian 
fluid, on purpose neglecting other complications present 
in yield stress fluids and thixotropic fluids, such as elas- 
tic and yield stresses and critical shear rates. Spatial 
resolution, necessary for the transient shear banding, is 
obtained when the rheological model is coupled to the 
Navier-Stokes equation for laminar flow in a concentric 
cylinder Couette device. The transient shear banding 
here occurs due to the shear rate - viscosity coupling in 
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the Navier-Stokes solution. We will demonstrate that 
such transients initiate since the shear stress exhibits a 
finite gradient in any rotational geometry. This gradient 
is amplified by the interaction with the shear thinning 
fluid. In the model, the timescales associated with the 
relaxation of the structure depend on the shear rate. A 
higher shear rate implies faster dynamics creating, with 
the coupled Navier-Stokes equations, an amplifying effect 
during the startup flow. 

Model. — Most microscopic descriptions of the struc- 
ture of a complex fluid are purely phenomenological, and 
based on a kinetic relaxation equation. The simplest A- 
models stipulate an evolution equation for a structural 
parameter, A j^iH]. This describes the internal order, 
such as the state of aggregation in colloids or the align- 
ment of particles with shear Usually, shear and /or 
temperature influences the temporal evolution of A [24| . 
The macroscopic rheology is obtained coupling A to a 
constitutive equation. Variants of A-models can be used 
to describe the flow curves of thixotropic, simple yield 
stress fluids, and shear thinning fluids [J, |25|. They gen- 
erally assume homogeneous flow, and thus are not appli- 
cable to shear bands. 

In aggregating suspensions and microgels a portion 
of the liquid is trapped due to the presence of solid 
structures [2^ 27 1. This is infact also the interpreta- 



tion shared in the literature for carbopol microstructure, 
which is observed to be formed of elastic sponge-like el- 
Therefore, instead of using an abstract 
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structure parameter (with no direct physical interpreta- 
tion) it is more sensible to select the immobilized volume 
fraction (p, which describes the jammed fluid, as it has 
been done in approaches using the Population Balance 
Equations (2^. The simplest time-dependent rheology 
follows by including a temporal relaxation of the volume 
fraction to a steady state. Unlike the relaxation rate, the 
steady state volume fraction is independent of the applied 
shear rate corresponding to a Newtonian fluid. Such a 
model can be thought as a Taylor-expansion of the dy- 
namical equation for the volume around the steady-state. 
The resulting kinetic equation is given 



dt {p^/p^oY 



B, 



(1) 



where As {Bg) is the kinetic constant for the shear growth 
(destruction) , 7 is the magnitude of the shear rate, k and 
7o (set to unity) both relate to the volume fraction sen- 
sitivity to shearing. A^, jio, and m describe the growth 
of jammed volume fraction due to the shear independent 
motion of the structure elements. The special case of 
m = k, presents a simple yield stress fluid, m > k gives 
a non-monotonic flow curve indicating thixotropy, and 
m < k produces shear thinning behavior. In what fol- 
lows, we flx the parameters Ag and Bs to 0.665 and 1.0, 
respectively. Since we are concentrating on a minimum 



model showing transient shear banding, in the follow- 
ing we set Ah = 0, and parametrize the initial volume 
fraction instead of specifying the sample history and the 
associated parameters. It is well known that in any prac- 
tical experiment the initial state depends on the sample 
history, and the shear independent structure dynamics. 
Therefore special attention has to be payed to the mea- 
suring protocol. The shear independent terms dominate 
the structure evolution at small shear rates making the 
flow curve non-Newtonian, but have negligible influence 
in the transient shear banding regime. There, it is rea- 
sonable to assume that both, the structure growth and 
destruction terms are determined by the shear rate jus- 
tifying such approximation. The exact form of the rate 
kernels are presently unclear and here we choose a sim- 
ple power-law dependence. The scaling factors 70 and /io 
make the equation unitless, and could be incorporated to 
the kernels Ag, Ai,, and Bg equally well. 

We use the well-known Krieger-Dougherty constitutive 
equation [29l|, that reads 



(2) 
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where (f>rn is the jamming volume fraction, /xq the liquid 
viscosity, and rj gives the scaling of the viscosity. From 
now on these quantities are set to = 0.68 (random 
sphere packing), /io = 1 mPas (water), and ry = 1.82 (3d| . 
These values relate to the initial conditions at which the 
transient shear banding appears, but are irrelevant e.g. 
for the scaling relations of the fluidization times 

The steady state volume fraction is </)ss 
consequence of leaving the shear independent terms out 
of the kinetic equation, the resulting steady state volume 
fraction is independent of the shear rate. As mentioned 
before, these shear independent processes are the ones 
that drive the system out of steady state implied by Eq. 
(1) when not sheared. This effect is incorporated to the 
model by simply initializing the system to the desired vol- 
ume fraction before starting the structure evolution. For 
shear thinning fluids having no yield stress, the present 
approximation is valid for all shear rates. Fluids hav- 
ing a yield stress, are properly described in the range 
where shear localization does not appear i. e. the flu- 
idization is complete. Besides flxing the steady state, Ag 
and Bg determine the rate of relaxation. Plugging <j)gg 
to the constitutive equation and applying the Newtonian 
assumption Ggg — ^{(pgg)^, gives the Newtonian steady 
state flow curve in Fig. [TJ The stress shows time and 
shear rate dependent exponential relaxation (the Inset), 
when started from (j)o 7^ 0ss , since the structure evolution 
follows Eq. H]). 

For the spatial resolution we couple the model with 
a continuum description of the fluid flow in a Couette 
rheometer geometry jsj]. The incompressible Navier- 
Stokes equation for laminar flow in that has the analytical 
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Figure 1. (Color online) The steady state flow curve, showing 
Newtonian response, with As — 0.665, Bs — 1.000. Inset: the 
transient behavior of the viscosity, ijio — 0.6799, k = 1.5. 



one-dimensional (radial) solution [32j 



7(r) 



dr 



(3) 



where r is the radial distance from the cell center, j{r) 
is the local shear rate, ila {^b) is the angular velocity of 
the inner (outer) cylinder, and Ra (Rb) is the radius of 
the inner (outer) cylinder. Eq. ([3]) implies that only the 
relative angular velocity of the two cylinders matters. In 
order to mimic the experiments of Ref. [12] only the inner 
cylinder is rotated and the radii are set to i?a = 23.9 mm 
and Rb = 25 mm. The Eq. ^ implies a spatial depen- 
dence for the volume fraction 0(i, r). To this end, a radial 
discretization is applied with a uniform grid with N sam- 
pling points, at each of which a separate cj) is evolved. 
Thus, the time evolution of the local flow field is ob- 
tained when the Equations ([T]) and @ are iteratively 
solved using a forward Euler algorithm under a constant 
global shear rate. Solving the structure and flow field 
evolutions this way assumes that the time-scales of the 
inertial effects and the structure kinetics is very different. 
The fiow field quickly adapts the changes caused by the 
slowly varying time-dependent viscosity. 

The global shear rate is defined as the radial average 
shear rate J^*" 7(r)/(i?f, — Ra) for a Newtonian fluid in 
the Couette geometry. This reads 



(7) = (rib - ^a) 



2RaRb 

Rb ^ ^1 



(4) 



Other definitions for the "engineering" shear rate are also 
used [3^ H^l- AH such are linearly proportional to the 
difference of the angular velocities, ilb — Oa. 

Results. — Varying the initial volume fraction and the 
kinetic exponent k uncovers three different startup flow 



scenarios as illustrated in the schematic phase diagram 
shown in Fig. [2j An almost homogeneous relaxation of 
the volume fraction profile arises at initial volume frac- 
tions above the steady state value and low kinetic expo- 
nents (Fig. [2Ia)). Increasing either makes the transient 
shear bands to appear (Fig.[21[b)). At initial volume frac- 
tions below the steady state, there is no transient shear 
bands (Fig. EJc)). The boundaries of that phase (grey 
area) are qualitative, as the transition from the homoge- 
neous relaxation Fig.[2l^a) to the transient shear banding 
Fig. [2Ib) occurs smoothly. 

t 
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Figure 2. (Color online) A schematic phase diagram of the 
instability /stability regimes as a function of k and (po- In 
regions (a) and (c) spatially homogeneous relaxation is ob- 
served, whereas region (b) exhibits transient shear banding. 

For a given value of k, the initial value 0o is chosen in 
order to be in the transient shear banding region (grey 
area in Fig. [5]). Thus, the simulation starts close to jam- 
ming, at (j)o = (j){r. t = Q) ^ Fig.[3l^b) shows the local 
shear rates corresponding to the evolution of the volume 
fraction plotted in Fig.[3l^a), from the numerical solution 
of Eqs. ^H^. The shear rate profiles show the develop- 
ment of a transient shear bands having two clearly dis- 
tinct bands evolving towards a homogeneous flow. Com- 
paring the velocity profiles, FiglUc), with the carpobol 
gel experiments 12| shows close similarities, even if the 
model here exhibits no yield stress. 

The origin of transient shear banding is a self-feeding 
mechanism, which can be understood considering the 
Eqs. ([l][3]). When the fluid starts to evolve at a homoge- 
neous high volume fraction, the 4>{r) decreases fastest at 
the regions, where j{r) has the largest value (Eq. (jlj): 
close to the inner cylinder of the device. The accelerated 
decrease of /i(r) (Eq. ^) due to the faster relaxation 
rate, further increases the j{r) (Eq. [3]) at the same loca- 
tion. This accelerates the decrease of the 0(r) (Eq. ([1])) 
at the same position creating a self- amplifying mecha- 
nism for the growth of transient shear banding. Since 
the steady state viscosity is constant, the (/)(r) will de- 
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Figure 3. (Color online) The volume fraction (upper panel) 
and corresponding shear rate (lower panel) during the startup 
flow. As = 0.665, Bs = 1.000, k = 1.5, (7) = 1.0 s"^ and 
(j>{r,t^ 0) = 0.6799. 



crease elsewhere in the device, only with slower rate due 
to lower j{r). Finally a homogeneous steady state pro- 
file is reached. If the feedback either from <j>{r) to j{r) 
(Eqs. (El) and or from j{r) to <j){r) (Eq. (P)) is 

not strong enough, the transient shear banding does not 
appear. The intensity of the feedback is adjusted by the 
exponent k and the derivative of the Eq. ([2]) at the cor- 
responding (j>{r). If (j>{r) at startup is below the steady 
state one, such a feedback loop does not exist; the growth 
of the (/)(r) is fastest at the high shear rates, thus pro- 
moting the increase of the /i(r) and the decrease of the 
7(r) at high shear rate regions. This reduces the shear 
rate differences in the gap. 

It has been observed experimentally, that in simple 
yield stress fluids the fluidization time decays with the 
global shear rate following a power-law. Similar be- 
haviour is observed here (Fig. 3]). Here, the power-law 
decay follows Tf ~ (7)"*^, with k fixed by/to the expo- 
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Figure 4. (Color online) The fluidization times against (7). 



nent k in Eq. ([T]). The other parameters or even the 
form of the viscosity function (Eq. ^) do not influence 
the power-law slope. The rest of the model parameters 
and the gap width simply change the vertical position of 
the resulting power-laws. 

The combination of Eq. ^ with = and Eq. ^ 
that describes the fluid in a Couette can be rewritten as 



70 
(7) 



dt 



iA,-B,4>ir,t))[F{r,t)f 



(5) 



where [F(r, t)] represents the geometry effects, and 
in particular is crucial for creating the transient shear 
bands. This form reveals a natural way to rescale the 
time: ^(7)'=. The measures dependent on the angular 
velocity can be rescaled with the global shear rate. 

Other ways to illustrate transient shear banding, be- 
sides the velocity profiles, are to plot the temporal evo- 
lution of the local shear rate at the gap edges, and the 
band width [l3|. These quantities are plotted with the 
same rescaling of time and shear rate in Fig. [5] (a) and 
(b), respectively. The band edge is estimated as the posi- 



tion, where 



— 0. The Fig.Oshows that there is 



a short induction period at small times, during which the 
shear rate localizes near the rotor. During the relaxation 
period, the shear rate decreases towards the steady state 
value as the transient shear bands vanish. 

Divoux et al. (2^ demonstrated that the fluidiza- 
tion time exponents in stress and shear rate controlled 
experiments relate to the steady state power-law expo- 
nent. Here, this equals unity. As one would expect, 
the fluidization for the stress controlled case follows a 
power-law tf ^ as does the shear rate controlled 
one tf A structural kinetics model producing 

Herschel-Bulkley flow curve with a control of the expo- 
nent can be speculated to change also the the fluidiza- 
tion exponents towards the experimental ones |22i] . To 
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Figure 5. (Color online) The normalized shear rate evolution 
close to the inner and outer cylinder (upper right panel), and 
the normalized evolution of the shear band edge position in 
the gap. 



the transient shear bands due to the shear dependent re- 
sponse of the fluid. Experimental observations of stress 
signatures associated to transient shear banding found 
in a cone-and-plate geometry, where the stress gradient 
is extremely small, also support this finding [36i]. Fur- 
thermore, in Eq. [3]we neglected the inertial terms in the 
Navier-Stokes equation, which could play a role in the 
start-up flow in such situations (referring to cone-and- 
plate). The self- amplifying mechanism between the flow 
and the fluid structure would be the same, but the origin 
of the shear inhomogeneity would come from the inertial 
terms rather than the flow geometry, which was the case 
in the circular Couette studied here. 

The kinetic exponent related to the dependence of (f) 
on 7, is connected to the fluidization exponent, in both 
the stress and shear controlled cases. This indicates that 
studying the fluidization experimentally gives detailed in- 
formation of the timescales of the internal relaxation pro- 
cesses in time-dependent fluids, which could be utilized to 
build proper structural models based on the experimental 
fluidization data. Such models could be established for 
instance around rheological models describing the volume 
fraction using Population Balances, which take into ac- 
count the particle size distribution and concentration as 
rheological parameters [26|. Presently, work is devoded 
along these lines of research. 
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the best of our knowledge, such a model does not ex- 
ist. Studying the properties of more advanced structural 
kinetics equations on the fluidization is therefore left to 
future studies. 

Conclusions. — We have studied transient shear band- 
ing, or fluidization of a simple time dependent fluid. We 
constructed a minimal model which lacks further com- 
plications such as the (visco)elasticity of the structure, 
a physically-motivated yield stress, the normal stresses 
suggested very recently |35j| to play a role in steady state 
shear banding [ill, or elastic stress overshoots Our 
analysis indicates that such conditions are not required 
for transient shear banding. On the contrary, transient 
shear banding should be a general feature of complex flu- 
ids like colloidal suspensions and microgels. 

The transient shear bands during the fluidization was 
found to originate from the initial shear inhomogeneity 
here arising from the Couette geometry, which is am- 
plified by the positive feedback coming from the shear 
dependent relaxation rates of the fluid. This mechanism 
should be present in all similar scenarios of fluidization: 
studying the couplings built into the model implies, that 
the same effect should occur in all practical measuring 
geometries for shear thinning complex fluids. This is 
since even the smallest stress gradient is enough to trigger 
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